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Abstract: We consider an inflationary model based only on renormalizable superpotential 
terms in which a superheavy scale F-term hybrid inflation (FHI) is followed by a Peccei-Quinn 
(PQ) phase transition. We show that the field which triggers the PQ phase transition influences 
drastically the inflationary dynamics and that the Universe undergoes a secondary phase of 
reheating after the PQ phase transition. Confronting FHI with the current observational data we 
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to avoid the gravitino problem. 
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1. Introduction 

In this talk, which is based on Ref. [|l|], we describe how we can achieve a cosmological 
scenario in which a superheavy F-term hybrid inflation (FHI) is followed by a Peccei-Quinn phase 
transition (PQPT) using two similar renormalizable superpotential terms. Below, we first briefly 
review the basic ingredients of our construction in Sec. [O] and Sec. L2 and outline the structure 



of our proposal in Sec. 1.3 



1.1 F-TERM Hybrid Inflation 

One of the most natural, popular and well-motivated inflationary model is the supersymmetric 
(SUSY) FHI [Q, ^, ^. It can be realized adopting the superpotential 

Wfhi = k-5(<1><I>-M2), (1.1) 

which is consistent with a continuous /^-symmetry [||] under which 

S e'"S, ^ OO, IVfhi ^ e'^Wmi- (1-2) 

Here, 5 is a left-handed (LH) superfield, singlet under a grand unified theory (GUT) gauge group 
G; <l> and <I> is a pair of LH superfields belonging to non-trivial conjugate representations of G, and 
reducing its rank by their vacuum expectation values (v.e.vs); K and M are parameters which can 
be made positive by field redefinitions. 

The SUSY potential induced by Wfhi in Eq. along the D-flat direction |<l>| = |<I>| is 

VPHI = fc^ \^-M^f + K^\S\^ (|<l>|2 + lop) . (1.3) 

W^FHi gives rise to FHI, since there is a F-flat direction, with <l> = = and constant potential 
energy Vp ~ K^M^, which is a local minimum of for S > M. Also, IVfhi leads to the spontaneous 
breaking of G, since the SUSY vacuum lies at 

(5)=0 and |(<1>)| = |(<I>)|=M, (1.4) 

with the non-zero v.e.vs of O and <I> developed along the Standard Model (SM) singlet directions. 

One of the shortcomings of FHI is the tension which, in general, exists between the pre- 
dicted (scalar) spectral index «s and the recent seven-year results ^ from the Wilkinson microwave 
anisotropy probe (WMAP7) satellite. Indeed, it is well-known that the realization of FHI within 
minimal Supergravity (SUGRA) leads to ?is which is just marginally consistent with the fitting of 
the WMAP7 data by the standard power-law cosmological model with cold dark matter and a 
cosmological constant (ACDM). One possible resolution of this problem is the addition to the 
Kahler potential of a non-minimal quatric term of the inflaton field with a convenient choice of 
its sign. As a consequence, a negative mass term for the inflaton is generated. In the largest part 
of the parameter space, the inflationary potential acquires a local maximum and minimum. Then, 
FHI of the hilltop type can occur as the inflaton rolls from this maximum down to smaller val- 
ues. Therefore, can become consistent with data, but only at the cost of an extra indispensable 
mild tuning of the initial conditions. Another possible complication is that the system may get 
trapped near the minimum of the inflationary potential and, consequently, no FHI takes place. 
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1.2 SUPERSYMMETRIZING THE PQ SOLUTION TO THE STRONG CP PROBLEM 

Due to the non-perturbative structure of the vacuum of SU (3)c the lagrangian of quantum 
chromodynamics (QCD) includes a CP- violating term, involving the strong coupling constant, g3, 
the gluon field-strength tensor, ^, and its dual, I.e., 

^QCD = 3II2 ^ + • • • With e<5- 10-'«, (1.5) 

since 6 is involved in the computation of the neutron electric dipole moment which is experimen- 
tally determined, with result 

^ 4.5 • 10"'^ < 2.9 • 10^2^ e-cm at 90% c.l. (1.6) 

The smallness of 6 consists the infamous strong CP problem. The most promising solution, pro- 
posed ^ by Peccei and Quinn, is to introduce a global color anomalous ?7(1)pq symmetry which 
is spontaneously broken at an energy scale fa ^ (10^ — 10'^) GeV, known as PQ energy scale. The 
Goldstone boson, a{x), associated with such symmetry breaking is called axion. The Lagrangian 
term resulting after the spontaneous symmetry breaking of the ?7(1)pq symmetry reads: 

^a = ^d^'ad^a + Ca^j%y%^,y, (1.7) 

where Ca is a model-dependent parameter. When considering the total lagrangian parts of Eqs. ( pTsj ) 
and (1.7), an effective potential for a appears, whose minimum is reached when the so-called 
(axion) misalignment angle vanishes, i.e., 

e = d + Ca^=0 or <a>=-d—- (1.8) 

Therefore, minimizing the potential with respect to a sets the offending CP- violating term to zero. 
Essentially, Q is promoted to a dynamical variable that evolves to its CP-conserving minimum, 
0=0, where Q can be seen as the phase of a new complex scalar field, named PQ field. 

Within a SUSY framework, the spontaneous breaking of ?7(1)pq can be obviously realized 
adopting a renormalizable superpotential, VKpq, similar to that of Eq. ( |1.1[ ) where S is replaced 
by another G and PQ singlet LH superfield, P, with the same R charge, while <I> and <l> are replaced 
by a pair of G singlet oppositely PQ-charged LH superfields, Q and Q. Indeed, the superpotential 

WpQ = JQ,P((2!2-/«V4), (1.9) 

is invariant under the ?7(1)pq transformations 

P ^ P,Q^ e'^Q^Q ^ e-'^Q (1.10) 

and lead to the F-term SUSY potential 

VpQF=K-2|(2e-/„V4|' + K-2|Pp(|e|2 + |2|2), (l.H) 

from where we can infer that U{1)pq can be spontaneously broken due to the following v.e.vs: 

(P) ~ and = /„ with 22 = 2(2 = 02 (1-12) 
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- since the sum of the arguments of {Q) and (Q) must be 0, Q and Q can be brought to the real axis 
by an appropriate PQ transformation. In reality, however, the total potential of the PQ fields is 

VPQ, = VPQF + K, where Va = ml{T)f^ {I - cos G) (1.13) 

(with Ma the temperature, T, dependent a mass) comes from nonperturbative QCD effects associ- 
ated with instantons [[To|], that break explicitly U (1)pq down to a discrete subgroup, where A'^ 
is the sum of the PQ charges of the SU (3)c triplets and anti triplets of the model. Therefore, the 



breakdown of by the v.e.vs in Eq. (1.12) may lead [10] to cosmologically catastrophic domain 



walls which, however, can be avoided 1 1 1 ] by introducing extra matter superfields - see Sec. 2.3 



A by-product of the ?7(1)pq spontaneous breaking is that we can achieve [12] a resolution 
of the /I -problem of MSSM by considering, e.g., a non-renormalizable superpotential term of the 
form Q^HiiH^/mp, which after the spontaneous breakdown of U (1)pq leads to the /i term of the 
MSSM, with ImI ~ I {Q)f /mp, which is of the right magnitude if | (2) | = fa/2 ~ 5 • 10" GeV 
and ~ (0.001 - 0.01) - here, mp ~ 2.44 • 10^^ GeV is the reduced Planck scale; //„ and Ha are 
the electroweak Higgses of MSSM which couple to up- and down-type quarks respectively. 

1.3 Outline 

The key point of our attempt in combining both ingredients (FHI and PQPT) described in 
Sec. [O] and L2 is that P can be regarded as the linear combination of the G x f/(l)pQ singlets with 



the 7? charge of the superpotential that does not couple to <txt> - cf. Ref. [|13|]. As a consequence, 
an unavoidable superpotential coupling SQQ and a.S — P mixing in the Kahler potential arise - see 
Sec. ^ These facts influence drastically the inflationary set-up described in Sec. ^. In addition, the 
value of P after FHI is to be kept larger than fa/2 so as to achieve an instantaneous domination 



of the PQ system over radiation in order to alleviate the gravitino (G) problem [14, 15]. These 
effects are presented in Sec. ^ We end up testing our model against observations in Sec. § and 
summarizing our results in Sec. ^. 

2. Model Description 

We below describe the structure of our model in Sec. \2.l\ we sketch its cosmological conse- 



quences in Sec. ^ and explain how we avoid the formation of domain walls in Sec. |23 



2,1 The General Set-up 

In order to explore our scenario, we identify G with the left-right symmetric gauge group 
Glr = SU (3)c X SU (2)l x SU (2)r x U(\)b-l, which can be broken down to the SM gauge group 
GsM = SU (3)c X SU (2)l x U (l)y through the v.e.vs acquired by a conjugate pair of SU (2)r dou- 
blet Higgs, <l> and <I>. As a consequence, no cosmic strings are produced in the end of FHI and, 
therefore, no extra restrictions on the parameters have to be imposed - cf. Ref. [[T^. The model 
possesses also three global U{\) symmetries. Namely, a (color) anomalous R symmetry U (1)r, an 
anomalous PQ symmetry U{\)^q and the baryon number symmetry U{\)g. The representations 
under Glr and the charges under the global symmetries of the various matter and Higgs superfields 
are presented in Table ^, which also contains n extra matter superfields (Da — and H^) required 
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2 


2 














//rf (1,2,-1/2) 











Table 1: The representations, the transformations under Glr, the decompositions under Gsm as well as 
the extra global charges of the superfields of our model. Here, Uq € SU{3)c, J/l € SU{2)i^, Ur £ SU (2)r 
and T and * stand for the transpose and the complex conjugate of a matrix respectively. 



for evading the domain-wall problem associated with PQPT together with a new imposed global 
U{l)jj symmetry - see Sec. |2.3[ 

In particular, the superpotential, W, of our model reads: 

W = Wfhi + WpQ + XSQQ + Wmssm + Wdw, (2.1) 
where Wfhi and Wpq are given by Eqs. (1.1) and ( |1.9| ) respectively and the anticipated in Sec. 1.3 
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unavoidable coupling is included. In addition, 

• VKmssm is the part of W which contains the usual terms of the Minimal SUSY SM (MSSM), 
supplemented by a mass term and Yukawa interactions for right-handed neutrinos, vf : 

W^MSSM = Am- \-yyij ^ +yiijlihl) +yqijqihq) . (2.2) 

Imp mp 

Here, the ith generation SU (2)l doublet LH quarks and leptons are denoted by qj and /,■ respectively, 
whereas the SU (2)r doublet antiquarks and antileptons by q'i and /[ respectively. The electro weak 
Higgs are contained in a SU (2)l x SU (2)r bidoublet Higgs h. The first term in the right-hand side 



(RHS) of Eq. generates the jj. term of MSSM via the PQ breaking scale - see Sec. [0|-, while 
the second term generates intermediate scale masses for vf and, thus, seesaw masses [||] for the 
light neutrinos - the coupling constant matrix yyij is considered diagonal. 



• Wdw is the part of W which gives intermediate scale masses via (Q) - see Sec. [ij - to 
Da — Da and Ha- Namely, 

IVdw = ^oMD^D, + Xh.QHI (2.3) 

where the coupling constant matrices Xoa and A/ja are considered diagonal. Although these matter 
fields acquire intermediate scale masses after the PQ breaking, the unification of the MSSM gauge 
coupling constants is not disrupted at one loop. In fact, if we estimate the contribution of D^^D^, 
and Ha to the coefficients b\, b2, and controlling ^VT\\ the one loop evolution of the three gauge 



coupling constants ^1,^2, and g3, we find that the quantities bj — b\ and Z73 — bi (which are |17] 
crucial for the unification of g\,g2^ and ^3) remain unaltered. 

The Kahler potential for our model can include interference terms of S and P even at the 
quadratic level, i.e., it has the form 

K = \S\^ + \P\^ + a{SP* +S*P)+b^-^ + c^-\ + d ^ ^ ' ' + ' ' ' {SP*+S*P) 

4mp 4mp mp 2mp 

+ •••, (2.4) 



+ 



4ml 



{SP*f + {S*Pf 



where all the coefficients a,b,c,d,e,f and g are taken, for simplicity, real. The ellipsis represents 
terms involving the waterfall fields (<I>, <l>, Q, and Q) which have negligible impact on our analysis. 

2.2 The Cosmological Scenario 

The F-term SUGRA scalar potential, Vsugra of our model can be found by applying the 
well-known formula - see e.g. Ref. [Q]: 

Vsugra = e"/'"'^ (f*, (Kji,)-' Fj - with F,- = Wj + K~- (2.5) 

Here, a subscript ,/ [, /*] denotes derivation with respect to the complex scalar field / [/*]. Taking 
the limit mp — )• 00, we can obtain the SUSY limit of Vsugra> Vf» which turns out to be 

Vf = j^^^(\K{^^-M^)+?iQQ\^+ k-2|(2!2-M2|')+ jc2|5|2(|<l>|2 + |<I,|2) 

+ \XSQ+KaPQ + ?iDaDaDa + XH,Ha\^+ |(2|'(|A5+JQ,Pp+A^a(|5a|' + |£>a|2)+A^al^a|' 



(l-a2 



T^a{Q*Q*-M^,) [k - M^) + XQQ]+ c.c. 



(2.6) 
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Super- 
fields 


SCALARS 
(12 REAL) 


Mass 
Squared 


Fermions 
(6 Weyl) 


Mass 
Squared 


<l>, O 

Q,Q 


Re&lm [V^±VS] 
Re&lm [e^±4] 
Re&lm [Q±Q] 


1^ |2 _i_ K{k-aK„)M^ 


Vv± - ^ 

Ve±Ve 


\Oa\' 



Table 2: The SUSY-breaking mass spectrum along the inflationary trajectory ofEq. (2.8a). Here, y/^ with 
x~v,V,e,e,Q and Q denote the chiral fermions associated with the superheldsx = V^, vj,e^,e^, 2 and Q 
respectively. 



where the complex scalar components of the superfields are denoted by the same symbol. From 
the potential in Eq. (2.6) and taking into account that fa,vje find that the SUSY vacuum Ues 



at the directions - cf. Eqs. (1.4) and (1.12): 



{S) ~ 0, (P) ~ 0, {Da) = {Da) = {Ha) = {^^) = {e%) = 0, 
(v|,) = (v|,)=M and (0e)=/«> 



(2.7a) 
(2.7b) 



where we have introduced the canonically normalized scalar field 0g = IQ = 2Q - cf. Eq. ( |1.12| ). 
As a consequence, W leads to a spontaneous breaking of Glr and U{1)pq. In addition, W gives 
rise to a stage of FHI and a PQPT, since Vp possesses two D- and F-flat directions for 



and S = e 



Q = Q = Da=Da=Ha=0 



e = e = and v|, = v| 



with a constant potential energy density respectively 

(a) Vmo ^ »c2mV(1 - a^) and (b) Vpqo = K^f^/16. 



(2.8a) 
(2.8b) 



(2.9) 



By constructing the scalar spectrum along the direction of Eq. ( 2.8a| ) - see Table I- we can deduce 
that it can be used as inflationary path since it corresponds to a classically flat valley of minima for 



\S\ > 



M 



and (b) \Oa\ > \l ^ — where Oa = XS+K,P- (2.10) 
1 — <i 



Since Vpqo <^ Vum, Vpqo can dominate over radiation after the end of FHI leading to a PQPT. This 
cosmological scenario can be attained if Eq. ( 2.1(]| a) is violated before Eq. (2. ICb), since, in this 
case, we obtain v|, = v|, = M and Q = Q = and not v|, = v|, = and Q = Q = fa/^. 

2.3 Evading the Domain-Wall Problem 

Soft SUSY breaking and instanton effects explicitly break U{1)rxU (1)pq to a discrete sub- 
group, which can be found, for every n, by solving the system of equations: 



^irR(W) ^ ^ 



4r = (mod 2n) 

-12r + 2(«-6)/7 = (mod 2n) 



where 



e^pPQ G [/(l)pQ, 



(2.11) 
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with r [p] being a U{1)r [?7(1)pq] rotation and the sum over / is appUed over all SU (3)c 3 and 3 of 
the model. We conclude that the unbroken subgroup is Z4 x Z2(„_6). It is then important to ensure 
that this subgroup is not spontaneously broken by (Q) and (Q), i.e., the equations 



e^'P'iQ) = {Q) and e-^'P-^{Q) = {Q) 2p, = (mod Iti) . 



(2.12) 



are satisfied identically - otherwise, cosmologically disastrous domain walls are produced [ |10| ] at 
PQPT. This goal can be accomplished by choosing n = 5 or n = 1. Therefore, for these n's, the 
domain-wall production during PQPT can be eluded. 

3. The Inflationary Era 



Below, we describe the saUent features of the inflationary potential in Sec. 3. 1 and we analyze 



the inflationary dynamics in Sec. |3.2. 
3.1 The Inflationary Potential 



The inflationary potential along the trajectory of Eq. ( [2.8a| ) can be written as 

Vm = Vmo + Vuis + Vnic , where 



(3.1) 



Vhio is the dominant contribution to Vhi along the F-flat direction, given in Eq. (2.9a) 



Vhis is the SUGRA corrections to Vhi which can be found by expanding Vsugra in Eq. (2.5) 

Vhio 



along the trajectory of Eq. (2.8a). Namely, 



Vms ^ 



:i-a2)m2 



Ai\S\^+An{S*P + PS*)+A2\P\ 



+ 



4(1 



Bi\S\^+B2\P\'^ 



+ Bi\S\^\P\^+{B4\S\'^+B5\P\^){S*P + PS*)+BJ{S*Pf + {P*S) 



Vhio / 2 / 2 , 2\ , 2 2\ , 



(3.2a) 



where the coefficients Ai — A3 and Bi —B^, given in Ref. [[T]], are functions of the coefficients a, ...,g 
in Eq. (2.4); the real fields o,s and q are the eigenvectors (corresponding to the eigenvalues m\) 



of the matrix involved in the quadratic part of Vhis- This can be worked out [|lj] after the quadratic 



part, Ksp, of K in Eq. ( |2.4t ) has been brought into a canonical form, i.e., we obtain also 

Ksp = \S\^ + \P\^ + a{SP* + S*P) = (a^ + / + q^) /2. 

Since ~ 0, its corresponding eigenvector, a, can be qualified as the inflaton. Note that we need 
the higher order terms of K in Eq. ( 2^ ) so that we obtain wPl < and therefore, observationally 
acceptable ?is's - see Sec. ^ Indeed, for a 7^ and b = c = d = e = f = g = we get = 0. 

• Vhic represents the contribution to Vm from one-loop radiative corrections, due to SUSY- 
breaking mass spectrum presented in Table 0, which can be calculated [|l^] to be 



Vmc ^ 



k"^Vhio 



^"(T3^ + 2j+ 16;r2(l-a2) r i^-^'W +2.' ' ^^"''^ 



87r2(i-a2) 

with X = |5|2(1 — a^) /M^ and Xa = \aa\^{l— a^) / k{X — aKa)M^. Here, we take into account that 
the dimensionality of the representations to which <l> and <I> {Q and Q\ belong is 2 [1] - see Table 1. 
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Figure 1: The evolution of a (black 
lines), s (gray lines), and q (light gray 
lines) as functions of N for the values of 
the parameters shown in Eq. (5.4), K = 
0.0023, a= -0.011,/? = -0.01 and /mi = 
3 • 10'^ GeV (solid lines) or fun = 2.5 • 
lO'^ GeV (dashed lines). For the employed 
values of parameters, the requirements of 
Sec. 5. 1 are fulfilled. 



3,2 The Inflationary Dynamics 

The equations of motion (e.o.m) of the various fields are ( ' = d/dt with t the cosmic time): 

f + 3Hf + Vuij = ^ H^f + 3//Y + Vuij = with f = o,s, and q (3.3) 

and ' = d/dN where N = ln (R/Rym). Here R{t) is the scale factor of the universe and the subscript 
"HIi" denotes values at the onset of FHI. We impose the following initial conditions (at N = 0): 



/hi; = /(O) = (1.5-4.5) • 10^^ GeV and /(O) = with f = a,s, or q. 



(3.4) 



When /he is large enough, s reaches an attractor and our results are independent of the precise value 
of /ho, as can be clearly deduced from Fig. |l|, where we plot a (black hnes), s (gray hnes), and 
q (light gray lines) as functions of N for fun = 3 • lO'^ GeV (solid lines) or /he = 2.5 • 10^' GeV 
(dashed lines). In both cases, we adopt the values of the parameters shown in Eq. (5^), K = 
0.0045, b = —0.01 and a = —0.011 which fulfill the requirements of Sec. 5A. For both choices of 
/hk's, we obtain ml = -0.0126, m^ = 1.83, n, = 0.968, A^hi* = 52, Omf = 4.08 • 10^^ GeV and 
i'Hif = 1.3 • 10^"^ GeV although in the first [second] case we obtain Nmi = 87.1 [A^hi = 175.5] -Nm 



and A'^Hi* are defined below Eq. (5.1) in Sec. 5.1. We observe that immediately after the onset of 
FHI, q decreases sharply, whereas the value of s at the end of FHI, suu, turns out to be just mildly. 



and not drastically reduced compared to (7 - in sharp contrast to the situation of Ref. [|19|]. This is 



due to the participation of s in both Eqs. ( |3.2a[ ) and ( |3.2q ). 



4. The Post-Inflationary Era 

We below describe the post-inflationary evolution of our model, presenting the dynamics of 
the two fields, a and s, in Sec. ^?l| and this of the two reheating processes in Sec. 4.2. For later 
convenience, we arrange in Table ^ the mass spectrum of our model at the SUSY vacuum of 



Eqs. (2.7a) and (2.7b) 



4.1 The Dynamics of Scalars 



When FHI is over, the inflaton system with mass mj - see Table ^ - consisting of the two 



complex scalar fields S and (5v|, — 5v^)/v2 - where = v, 



M and 5v^ 



-M - settles 
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Table 3: The mass spectrum of the model at the SUSY vacuum ofEq. (2.7a) and (2.7b). Here, A'^ [A^] 
with m = 1,2,3 are gauge bosons [gauginos] associated with the SU{2)^, while these corresponding to 
U{1)b-l are denoted byAB-i [Ab-l]- Also, y/^ with x = S, V, V,e,e,P, g, Q,Da,Da and i/a denote the Weyl 
spinors associated with the superBeldsx ^ S,V^,V^,e^^,e^^,P,Q,Q,Da,Da and H^, respectively. 



into a phase of damped oscillations and decays reheating the universe to a temperature 



/ 72 \ '^"^ 1 

^irh = ( ^ \/rimp, where H = -— A^mj (4.1) 



is the decay width emerging from the third term in the RHS of Eq. ( [2.1[ ). Here, gi^h* — 438.75 [girh* — 
5 13.75] for n = 5 [« = 7] counts the relativistic degrees of freedom of the model. 

For A ~ (0.05 — 0.1), we get Tirh > Vpgg. Therefore, we obtain matter domination (MD) for 



T > Tii-h and radiation domination (RD) for VpQQ < T < ^irh- During MD, s [|I9|,|2^ acquires an 
effective mass equal to y^3/2H. Solving its e.o.m for N > Nui, we can extract its value, ^'pgi, - 
and the corresponding value of P, PpQi - at T = Tirh which coincides with its value at the onset of 
PQPT since, during the subsequent RD era, s remains [11^, 20 1 frozen. Namely we find 



^Qi = ^P'^'PQi with SpQi 



Plrh 

Vhio 



1/4 



71 

Suit and pirh = ^q^ia*?!*! 



(4.2) 



where Ap is a function of the coefficients of K in Eq. (2.4). 



For T < VpQQ, W in Eq. ( [2. 1[ ) is dominated by Wpq in Eq. ( |1.9| ) and the relevant F-term scalar 



potential is given in Eq. ( |1.11[ ) which along the flat direction of Eq. ( 2.8b| ) gives rise to the con 



stant potential energy density of Eq. (2.9b). Assuming gravity mediated soft SUSY breaking, the 



potential along the direction of Eq. ( [2.8b| ) for |P| > fa/2 has the form 



VPQ VpQo + ml |Pp - 72^ |ap| \P\ + (in ^ + ^) , (4.3) 

where the 2nd and 3rd contributions arise from soft SUSY breaking effects and the forth contribu- 
tion represents the 1-loop corrections [^] due to the SUSY breaking Mainly due to this last 
contribution, Vpq does not give rise to another FHI, since the rj -criterion is spoiled. Nonetheless, 
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when l^l < fa/2, an instability occurs along the iPj-axis triggering thereby a PQPT. If, in addition, 
I^PQil > /a/2 we obtain an out-of-equilibrium decay of the PQ system, i.e., a secondary reheating. 

During this latter phase, the PQ system with mass mpQ - see Table |3| - comprised of the 
complex fields P and (52 + 52)/\/2- where 5Q = Q — fa/2 and 5Q = Q — fa/2 - enters a phase 
of oscillations reheating the universe to the temperature 



f 12 



2rh* 



1/4 



Y^r2mp, where r2 



fa 

2m\ 



(4.4) 



is the decay width emerging from the first term in the RHS of Eq. (2.2). Also, g2rh* = 232.5 counts 



the relativistic degrees of freedom of MSSM plus the content of the axion supermultiplet. 

4.2 The Dynamics of Reheating Processes 

A more accurate description of the reheating dynamics can be obtained by solving the relevant 
Boltzmann equations. In particular, the energy density, pi [pz], of the oscillatory system which 
reheats the universe at the temperature Ti^h [72rh], the energy density of produced radiation, Pr, 
and the number density of G, n^, satisfy the equations 



pi+3//pi+ripi =0, 

P2 + 3//p2 + r2P2 = 0, 

PR + 4HpR - Tipi - r2P2 = 0, 

0, 



with 



H = 



= 37ij:Ucigjln{k/gi) /I6i;{3)ml, 
= i;{^)Tyn\T = 30pR/g,n\ 



(4.5) 



Here, {ki) = (1.634,1.312,1.271), (c,) = (33/5,27,72) and g,(r) = gi** [g*{T) = g2A*] for T > 
TpQ [T < Tpq] where Trq is defined as the solution of the equation pR(rpQ) = Vpqq. We use the 
following initial conditions - the quantities below are considered as functions of the independent 
variable A^^ = In (/?//?Hif) with /?Hif being the value of the scale factor at the end of FHI: 



Pi (0) = Vhio, Pr(0) = nMO) = 0, and P2(A^pq) = Vpqo, 



(4.6) 



where A^^pq is the value of N corresponding to the temperature Tpq. 

In Fig. ^, we illustrate the cosmological evolution of the quantities log p, with / = 1 (dotted 
gray line), / = 2 (dashed gray line), and / = R (gray line), logVpQo (black dashed line), and logFg 
(black solid line) as functions of log T for the values of the parameters adopted in Fig. |TJ We 
observe that FHI is followed successively by a MD era, which lasts until T = 6 • lO^-' GeV ~ Tirh 
(where pi = Pr), a RD epoch, terminated at Tpq = 1.4 • lO^*' GeV, a MD era, completed at T = 
3.1-10'* GeV ~ r2rh (where p2 = Pr) and followed by the conventional RD epoch. We also see 
that the G abundance immediately after FHI is Y^^ = 5.5 ■ 10^^ which can be estimated by [14, 15] 



hh) 



1.9- 10^12 



-Irh 







with s : 



45 ^ 



(4.7) 



IQio GeV 

However, the G abundance decreases sharply to Y^^ = 1.6 • 10^*"* which can be approximated by 



30' 



Irh* 



1/4 



72rh 



Y ~ 

1/4 IG- 
PQO 



(4.8) 
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Figure 2: The evolution of the quan- 
tities logp, with i — 1 (gray dotted line), 
i = 2 (gray dashed line), i — R (gray line), 
logVpQo (black dashed line), and logFg 
(black solid line) as functions of \ogT 
for K = 0.0023, a = -0.011, b = -0.01 
and the values of the remaining parameters 
shown in Eq. (5.4). For the employed vales 
of parameters, the requirements of Sec. 5.1 
are fulfilled. 



log(T/GeV) 



1 /4 

We observe that Y^^ is suppressed relative to Y^^ by the ratio r2rh/VpQo ^ 1 due to the entropy 
released during the out-of-equilibrium decay of the PQ system. Interestingly enough, the dilution 
of is independent of jq, - see Eqs. ( |2.9| b) and (4.4). 



5. Testing Against Observations 



We below exhibit the constraints that we impose on our cosmological set-up in Sec. 5. 1 and 



delineate the allowed parameter space of our model in Sec. 5.2. 

5.1 Observational Constraints 

The parameters of our model can be restricted imposing the following requirements - note that 
in the point (v) below we adopt an updated, compared to our analysis in Ref. version of the 
relevant constraint : 



(i) The violation of the instability conditions in Eq. ( 2.1(\ ) occurs according to the desired order. 



(ii) The number of e-foldings A^hi* that the scale = 0.002/Mpc suffered during FHI has to be 
sufficient to resolve the horizon and flatness problems of Standard Big Bang cosmology: 



1/4 



1 



1/4 



A7 A7 A7 n , -1 HIO ^1 'PQO , ^1 ^lrh72rh 

Niiu = Nm —N^ ~ 23 + - In — — In — h - In ^ , 

3 1 GeV 3 1 GeV 3 1 GeV^ ' 



PQO 



(5.1) 



where A^^^ and N^i are the values of N from the onset of FHI until crossed outside the 
horizon of FHI and the end of FHI, respectively. Nui is the largest N at which we obtain 
violation of Eq. ( 2.10| a) or of the condition: 



max{e(a(A^)),|77(a(A^))|} < 1, with 



Vm,c 
Vm 



. 2 ^HI.fJCT 

and Tj ~ mp 



(5.2) 



(iii) The power spectrum of the curvature perturbation at k = k^ is to be confronted with the 
WMAP7 data: 



3/2 



HI 



2V3Kml\Vm.c 



4.93-10 



-5 



(5.3) 
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(iv) The mass, gM, of the Ughtest gauge boson at the SUSY vacuum - see Table 3 - is to take the 
value dictated by the unification of the gauge coupling constants within MSSM, i.e., 



~ 2 • 10^^ GeV ^ M ~ 2.86 • 10^^ GeV with g ~ 0.7, 



(5.4) 



being the value of the unified gauge coupling constant - not to be confused with the coeffi- 



cient g appearing in Eq. (2.4). Note that Glr is considered embedded in the 50(10). 



(v) The spectral index, n^, is to be consistent with the fitting of the WMAP7 results by the ACDM 
model (with negligible running — 0), i.e., 

n, = l-6e{N^) + 27] (A^,) = 0.968 ± 0.024 ^ 0.944 < < 0.992 at95%c.l. (5.5) 

(vi) In order for the PQPT to take place after a short temporary domination of VpQo, we require: 

I^PQil >fa/2 SpQi>fa/Ap. (5.6) 

(vii) Assuming unstable G, we impose an upper bound on order to avoid problems with the 



standard Big Bang nucleosynthesis [15|: 

10-14 

10-13 



Y ~< 



for G mass m; 



0.69 TeV 
10.6 TeV. 



(5.7) 



5.2 Numerical Results 

As can be seen from the analysis above, our cosmological set-up depends on the following 
parameters: k:, jq,, A, /o, A^, n, a, b, c, d, e, f, and g. We fix throughout our computation: 



A = 0.1, /„ = 10^2 GeV, Ka = Xa= 0.01, n = 5 and c = d 



:/ = g=0.1. 



(5.8) 



The chosen fa and result to ~ 1 TeV via the first term of the RHS of Eq. (^^. Also, the 
selected jq, and A play a crucial role in the determination of Ti^h and Tjrh - via Eq. (4J_) and ( pf ) 
and facilitate the violation of the conditions in Eq. ( 2.10D in the desired order. Their variation, 
thought, does not cause drastic changes in the inflationary predictions. The same is also valid for 



the fixed in Eq. (5.8) parameters of K, in Eq. (2.4) which - contrary to a and b - do not influence 
the computation of m\_ and . As we show below, the selected values above give us a wide and 
natural allowed region of the remaining fundamental inflationary parameters (k", a, and b). 

Besides the parameters above, in our computation, we use as input parameters the quantities N., 
and /he with f=a,s, and q. We set /he ~ (1.5 - 3.5) • lO" GeV so as to obtain A^hi '^10- 140. 
We then restrict M and A^* so that Eqs. ( |5.1[ ) and ( ^3| ) are fulfilled. It is gratifying that our model 



supports solutions which simultaneously fulfill Eqs. ( p.3[ ) and ( |5.4[ ) contrary to most realizations of 
cf. Ref. I 



FHI 



which requires, via Eq. (5.3), M's lower than those indicated in Eq. (5.4). We 



finally check if the preferred hierarchy in the violation of Eqs. (2. ICa) and (2.10 b) is achieved and 



proceed imposing the requirements (v) - (vii) of Sec. 5.1 



Letting a vary for a number of fixed values of b, we can depict the values allowed by all the 
constraints of Sec. ^ in the K — a plane - see the left plot of Fig. |[ The various lines terminate 
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b= 0.05 








b= 0.03 








b= 0.01 








b= - 0.01 
b = - 0.05 






1.5 


2.0 


2.5 3.0 

k(10-^) 


3.5 




0.006, K = 0.0037 
0.011, K = 0.0023 
0.0165, K = 0.00145' 



25 



30 



a (10'" GeV) 



Figure 3: Allowed values by the requirements of Sec. 5.1 in the K — a plane for various b's indicated 
on the curves (left) and the variation of Vhi as a function of a for b = —0.01 and {a,K)'s indicated in the 
graph (right). In both graphs we use the values of the parameters shown in Eq. (5.4). For the right graph 
we set /mi = 3 • lO" GeV (ns = 0.992, dashed line) or fun = 2.5 • lO" GeV («s = 0.968, solid line) or 
/hu = 2 • lO'^ GeV (hs = 0.944, dotted line). The values corresponding to a* and Of are also depicted. 



at low [high] k:'s due to the saturation of Eq. (5^) from below [above]. We readily conclude that 
the allowed {a,b)'s for fixed ?is are almost JC-independent. This is because ml is fixed too. In 
particular, for n, = 0.944, 0.968 and 0.992, we have -ml ~ 0.0179, 0.0125 and 0.0078 and K = 
0.00125,0.002 and 0.0037, respectively. In all cases, ~ 1.82, ~ 5 • 10"'' and ~ 2 • lO"^'*. 
Therefore, our scenario can be realized for both signs of a and b, contrary to the cases studied in 
Ref. [^ where negative b's are necessitated. Also, compared the extracted Y^^'s with the bounds 
of Eq. ( |5.7[ ), we infer that G with masses even lower than 10 TeV become observationally safe. 

One of the outstanding features of our proposal is that the reduction of can be attained 
without disturbing the monotonicity of the potential - cf. Ref. [^]. This fact is highlighted in the 
right plot of Fig. ^, where we present the variation of the inflationary potential Vhi as a function 
of a, for b = —0.01 and three pairs of a and k:'s, shown in the graph, corresponding to rig = 0.944 
(dotted line), 0.968 (solid line) and 0.992 (dashed line). The values corresponding to a* and Of are 
also designed. We observe that for large a's, Vhi develops an oscillatory behavior due to the initial 
oscillations of s and q - see Fig. |l|. However, Vui for lower a's remains monotonic and, therefore, 
no complications arise in the realization of FHI. 



6. Conclusions 



We showed that, combining FHI with a PQPT based on renormalizable superpotential terms, 
we can obtain: (i) Observationally viable FHI at the SUSY GUT scale with natural values, ±(0.01 — 
0.1), for the model parameters; (ii) a simultaneous resolution of the strong CP and /i problems of 
MSSM; (iii) a second stage of reheating after PQPT, which leads to observationally safe values of the 
G abundance. An important prerequisite for all these is that the field, which triggers PQPT, remains 
after FHI well above the PQ scale thanks to (i) its participation in the SUGRA and logarithmic 
corrections during FHI and (ii) the high reheat temperature after the same period. A noteworthy 
open issue of our scenario is this of baryogenesis which cannot be processed via non-thermal 
leptogenesis [Q] since the produced lepton asymmetry after FHI is efficiently diluted. 
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